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ABSTRACT 

In this Chapter we deal with the attempts to measure the general relativistic 
gravitomagnetic Lense-Thirring effect with the Satellite Laser Ranging (SLR) 
technique applied to the existing LAGEOS and LAGEOS II terrestrial satellites 
and to the recently approved LARES. According to general relativity, the rota- 
tion of the central body which acts as source of the gravitational field changes the 
spatial orientation of the orbit of test particles determined by the longitude of the 
ascending node Q. Such a shift, which is cumulative in time, is of the order of 2 
m yr~^ in the case of the LAGEOS satellites. Extracting this signature from the 
data is a demanding task because of many classical orbital perturbations having 
the same pattern as the gravitomagnetic one, like those induced by the centrifu- 
gal oblateness of the Earth which represents a major source of systematic error. 
The first issue addressed here is: are the so far published evaluations of the sys- 
tematic bias induced by the uncertainty in the even zonal harmonic coefficients 
Jg of the multipolar expansion of the Earth's geopotential reliable and realis- 
tic? The answer is negative. Indeed, if the difference AJ^ among the even zon- 
als estimated in different global solutions (EIGEN-GRACE02S, EIGEN-CG03C, 
GGM02S, GGM03S, ITG-Grace02s, ITG-Grace03s, JEM01-RL03B, EGM2008) 
is assumed for the uncertainties 6J£ instead of using their more or less cali- 
brated covariance sigmas aj^, it turns out that the systematic error 6fi in the 
Lense-Thirring test with the nodes of LAGEOS and LAGEOS II is about 3 to 
4 times larger than in the evaluations so far published based on the use of the 
sigmas of one model at a time separately, amounting up to 37% for the pair 
EIGEN-GRACE02S/ITG-Grace03s. The comparison among the other recent 
GRACE-based models yields bias as large as about 25 — 30%. The major dis- 
crepancies still occur for J4, Jg and Jg, which are just the zonals the combined 
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LAGEOS/LAGOES II nodes are most sensitive to. The second issue is the pos- 
sibility of reaching a realistic total accuracy of 1% with LAGEOS, LAGEOS II 
and LARES, which will be launched in the near future. While LAGEOS and 
LAGEOS II fly at altitudes of about 6000 km, LARES will be likely placed at 
an altitude of 1200 km. Thus, it will be sensitive to much more even zonals than 
LAGEOS and LAGEOS II. Their corrupting impact has been evaluated with the 
standard Kaula's approach up to degree £ = 60 by using AJ^; it turns out that 
it may be as large as some tens percent. The different orbit of LARES may also 
have some consequences on the non-gravitational orbital perturbations affecting 
it which might further degrade the obtainable accuracy. 



Introduction 



In the weak-field and slow motion approximation, valid when the magnitude of 
the gravitational potentials U and velocities v characteristic of the problem under 
examination are smaller with respect to the speed of light c, i.e. for U/c^,v/c <^ 1, the 
Einstein field equations of general relativity get linearized resembling to the Maxwellian 
equations of electromagntism. As a consequence, a gravitomagnetic field, induced by 
the off-diagonal components goi,^ = 1,2,3 of the space-time metric tensor related to the 



mass-energy currents of the source of the gravitational field, arises (IMashhoon 



2003); it 



has no classical counterparts in Newtonian mechanics. The gravitomagnetic field affects 
orbiting test particles, precessing gyroscopes, moving clocks and at oms and propagating 



electromagnetic waves (IRuggiero and Tartaglia 



20021 : 



S chafer 



2004 ) . Perhaps , the most 



famous gravitomagnetic effects are the precession of the axis of a gyroscope (jPugh 



1959 
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Schifj Il960l ) and the Lense-Thirring^ precessions (ILense and Thirring] Il918l ) of the orbit 
of a test particle, both occurring in the field of a central slowly rotating mass like, e.g., 
our planet. Direct, undisputable measurements of such fundamental predictions of general 
relativity are not yet available. 



The measurement of the gyroscope precession in the Ea rth's gravitat i onal f ield h as 



been the goal of the dedicated space-baseqj GP-B mission f lEveritt et al 



1974 



20011 ) 



launched in 2004 and carrying onboard four superconducting gyroscopes; its data analysis 
is still ongoing. The target accuracy was originally 1%, but it is still unclear if the GP-B 
team will succeed in reaching such a goal because of some unmodelled effects affecting the 
gyroscopes: 1) a time variation in the polhode motion of the gyroscopes and 2) very large 
classical misalignment torques on the gyroscopes. 

In this paper we will focus on the measurement of the Lense- Thirring effect in the 
gravitational field of the Earth. It consists of a secular rate of the longitude of the ascending 
node Q 



2GJ 



LT 



and of the argument of pericentre u 



c2a3(l-e2)3/2' 



GGJcosi 



(2) 



c2a3(l -e2)3/2' 

of the orbit of a test particle. In eq. ([1]) and eq. ([2]) G is the Newtonian constant of 
gravitation, J is the proper angular momentum of the central body, a and e are the 
semimajor axis and the eccentricity, respectively, of the test particle's orbit and i is its 
inclination to the central body's equator. The semimajor axis a determines the size of the 



^According to an interesting historical analysis recently performed in (jPfister 
would be more correct to speak about an Einstein-Thirring-Lense effect. 



20071), it 



See on the WEB http://einstein.stanford.edu/ 
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ellipse, while its shape is controlled by the eccentricity < e < 1; an orbit with e = is 
a circle. The angles Q and uj fix the orientation of the orbit in the inertial space and in 
the orbital plane, respectively. Q, u and i can be viewed as the three Euler angles which 
determine the orientation of a rigid body with respect to an inertial frame. In Figure [1] we 
illustrate the geometry of a Keplerian orbit. 

In this Chapter we will critically discuss the following two topics 



The realistic evaluation in Section [2] of the total accuracy in the test performed in 
rec ent years with the existing Earth's artificial sat ellites LAGEOS and LAGEOS 



II ( ICiufolini and Pavlis 



2004 



Ciufolini et al 



20061). LAGEOS was put into orbit 



in 1976, followed by its twin LAGEOS II in 1992; they are passive, spherical 
spacecraft entirely covered by retrorefiectors which allow for their accurate 
tracking through laser pulses sent from Earth-based ground stations according to 
the Satellite Laser Ranging (SLR) technique. They orbit at altitudes of about 
6000 km (olageos = 12270 km, olageos ii = 12163 km) in nearly circular paths 
(clageos = 0.0045, clageos ii = 0.014) inclined by 110 deg and 52.65 deg, respectively, 
to the Earth's equator. The Lense-Thirring effect for their nodes amounts to about 
30 milliarcseconds per year (mas yr~^) which correspond to about 2 m yr^^ at the 
LAGEOS altitudes. 

The idea of measuring the Lense-Thirring node rate with the just launched LAGEOS 
satellite, along with the oth er SLR targets orbiting at that time, was put forth by 



Cugusi and Proverbid (119781 ). Tests have starte d to be effectively pe rformed later by 



using the L AGEOS and L AGEOS II satellites (jCiufolini et al 



a strategy (jCiufolini 



19961 ). according to 



19961 ) involving the use of a suitable linear combination of the 
nodes Vt of both satellites and the perigee uj of LAGEOS II. This was done to reduce 
the impact of the most relevant source of systematic bias, i.e. the mismodelling in 



Fig. 1. — Keplerian orbit. The longitude of the ascending node is counted from a reference 
X direction in the reference {XY} plane, chosen coincident with the equatorial plane of the 
rotating body of mass M and proper angular momentum J, to the line of the nodes, i.e. the 
intersection of the orbital plane with the reference plane. The argument of pericentre uo is an 
angle in the orbital plane counted from the line of the nodes; the location of the pericentre is 
marked with 11. The time-dependent position of the moving test particle of mass m is given 
by the true anomaly /, counted anticlockwise from the pericentre; is an azimuthal angle 
in the {XY} plane. L is the orbital angular momentum, perpendicular to the orbital plane. 
The inclination between the orbital plane and the equatorial plane oi M is i. Courtesy by 
H. Lichtenegger, OEAW, Graz. 
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the even = 2,4,6...) zonal [m = 0) harmonic coefficients Ji of the multipolar 
expansion of the Newtonian part of the terrestrial gravitational potential due to 
the diurnal rotation (they induce secular precessions on the node and perigee of a 
terrestrial satellite much larger than the gravitomagnetic ones. The Ji coefficients 
cannot be theoretically computed but must be measured by using dedicated satellites; 
see Section [2]): the three-elem ents combination used allowed for removing the 



uncertainties in J2 and Jd.. In 



using thcl EGM96 flLemoine et al 



Ciu: 



blini et al 



1998a) a ^ 20% test was reported by 



19981 ) Earth gravity model; subsequent detailed 
analyses showed that such an evaluation of the total error budget was overly optimistic 

e computation of the total bias due to the even zonals 



in view of t 


le likely unreliab 


e 


( lorio 


2003; 


Ries et al 


2003a 


b) 



b|). An analogous, huge underest imation turned ou t to 



19871 ) hke 



hold also for the effect of the non-gravitational perturbations (iMilani et al 
the direct solar radiation pressure, the Earth's albedo, various subtle thermal effects 
depen ding on the physical properties of the satelhtes' surfaces and their rotational 



state ( 



Inversi and Vespe 



Lucchesi et al 



2004 



1994: 



Vespe 



Ries et al 



999 



Lucchesi 



2001 



Lucchesi 



2002 



2003 



2004 



2003al ). which the perigees of LAGEOS-like satellites 



are particul arly sensitive to. As a result, the realistic total error budget in the test 



reported in (jCiufolini et al 



1998bl ) might be as large as 60 — 90% or even more (by 



considering EGM96 only). 



^Contrary to the subsequent models based on the dedicated satellites CHAMP 



(http://www-app2.gfz-potsdam.de/pbl/op/champ/index_CHAMP.html) and GRACE 



(http://www-app2.gfz-potsdam.de/pbl/op/grace/index_GRACE.html), EGM96 relies 
upon multidecadal tracking of SLR data of a constellation of geodetic satellites including 
LAGEOS and LAGEOS II as well; thus the possibility of a sort of a — priori 'imprinting' of 
the Lense-Thirring effect itself, not solved-for in EGM96, cannot be neglected. 
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The observable used in (Ciufolini and 



GRACE02S model flReigber et al 



avlisl 120041 ) with the GRACE-only EIGEN- 



20051 ) was the following li near combination 



nodes of LAGEOS and LAGEOS II, explicitly com puted in (jlorio and Morea 



following the approach put forth in (jCiufolini 



oft 



he 



19961) 



20041 ) 



where 



Cl 



^LAGEOS 
^LAGEOS II 



/ = f^LAGEOS ^ cif^L^G'^OS 11^ 



/ 1 2 \ ^ / 

cos ^LAGEOS / J- — ClagEOS II \ / '^LAGEOS II 



(3) 



cos ZlAGEOS II 



1 - p2 

'^LAGEOS 



V C^LAGEOS 



7/2 



(4) 



19661) 



The coefficients ^l,e of the aliasing classical node precessions (iKaulal 
^Hass = ^j J e induced by the even zonals have been analytically worked out in, 



e.g. (jlorio 



20031 ): they yield Ci = 0.544. The Lense-Thirring signature of eq. 



amounts to 47.8 mas yr^^. The combination eq. ([3]) allows, by construction, to 
remove the aliasing effects due to the static and time-varying parts of the first even 
zonal J2. The nominal (i.e. computed with the estimated values of Ji, i = 4,6...) 
bias due to the remaining higher degree even zonals would amount to about 10^ mas 
yr~^; the need of a careful and reliable modeling of such an important source of 
systematic bias is, thus, quite apparent. Conversely, the nodes of the LAGEOS-type 
spacecraft are affected by the non-gravitational accelera t ions at a ~ 1% level of 



the 



jense-Thirring effect (ILucchesi 



2001 



Lucchesi 



2002 



2003 



2004 



Lucchesi et al 



20041 ). For a comprehensive, up-to-date overview of the num erous and su btle issues 



concerning the measurement of the Lense-Thirring effect see (jlorio 



2007df ) 



The possibility that the LARES mission, recently approved by the Italian Space 
Agency (ASI), will be able to measure the Lense-Thirring node precession with an 
accuracy of the order of 1% (Section [3|). 



^See also (jPavlis 



2002 



Ries et al 



2003 
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In (jVan Patten and Everitt 



1976a 



van Patten and Everitt 



1976bl ) it was proposed to 



measure the Lense-Thirring precession of the nodes of a pair of counter-orbiting 
spacecraft to be launched in terrestrial polar orbits and endowed with drag-free 



apparatus. A somew 



in 1986 by (jCiufolini 



lat eq uivalent, cheaper version of such an idea was put forth 



19861 ) who proposed to launch a passive, geodetic satellite in 
an orbit identical to that of LAGEOS apart from the orbital planes which should 
have been displ aced by 1 8 deg apart. The measurable quantity was, in the case of 



the proposal by ICiufolinil (jl986l ). the sum of the nodes of LAGEOS and of the new 
spacecraft, later named LAGEOS III, LARES, WEBER-SAT, in order to cancel to a 
high level of accuracy the corrupting effect of the multipoles of the Newtonian part of 
the terrestrial gravitational potential which represent the major source of systemat ic 



error 



989 



see Section [21) . Altho ugh extensively studied by various groups (IRies et al 



Ciufolini et al 



(llorio et al 



1998al ). such an idea was not implemented for many years. In 



2OO2I ) it was proposed to include also the d ata from LA GEOS II by using 



2005al ) to be potentially 



a different observable. Such an approach was proven in (llorid 
useful in making the constraints on the orbital configuration of the new SLR satellite 
less stringent than it was originally required in view of the recent improvements in 
our knowledge of the classical part of the terrestrial gravitational potential due to the 
dedicated CHAMP and, especially, GRACE missions. 

Since reaching high altitudes and minimizing the unavoidable orbital injection errors 
is expensive, it was explored the possibility of discardi ng LAGEOS and LAGEOS 



II using a low- altitude, nearly pola r orbit for LARES ( ILucchesi and Paolozzi 



2001 



Ciufolini 



20061), but in fllorio 



2002 



2007al ) it was proven that such alternative 



approaches are not feasible. It was also sugges ted that LARES would be able to probe 



alternative the ories of grayity ( 



be impossible ( llorio 



2005b 



Ciufolini 



2007bh . 



20041 ). but also in this case it turned out to 
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The stalemate came to an end when ASI recently made the following official an- 



nouncement ( http: / / www.asi.it / SiteEN/MotorSearchFullText.aspx?keyw=LARES ) : 



"On February 8, the ASI board approved funding for the LARES mission, that 
will be launched with VEGAs maiden flight before the end of 2008. LARES is a 
passive satellite with laser mirrors, and will be used to measure the Lense-Thirring 
effect." The Italian version of the announcement yields some more information 
specifying that LARES, designed in collaboration with National Institute of Nuclear 
Physics (INFN), is currently under construction by Carlo Gavazzi Space SpA; 
its Principal Investigator (PI) is I. Ciufolini and its scientific goal is to measure 
at a 1% level the Lense-Thirring effect in the gravitational field of the Earth. 
Concerning the orbital configuration of LARES, the ASI website says about VEGA 



that (http://www.asi.it/SiteEN/ContentSite.aspx?Area=Accesso-|-allo-|-spazio ): "[... 



VEGA can place a 15.000 kg satellite on a low polar orbit, 700 km from the Earth. 
By lowering the orbit inclination it can launch heavier payloads, whereas diminishing 
the payload mass it can achieve higher o r bits. [...]" In the latest communication to 



INFN, Rome, 30 January 2008, 



Ciufolinil (120081 ) writes that LARES will be launched 



with a semimajor axis of approximately 7600 km and an inclination between 60 
and 80 deg. More precise information can be retrieved in Section 5.1, pag 9 of 
the document Educational Payload on the Vega Maiden Flight Call For CubeSat 
Proposals, European Space Agency, Issue 1 11 February 2008, downloadable at 



http://esamultimedia.esa.int/docs/LEX-EC/CubeSat%20CFP%20issue%201.pdf. It 



is written there that LARES will be launched into a circular orbit with altitude 
h = 1200 km, corresponding to a semimajor axis olares = 7578 km, and inclination 
i = 71 deg to the Earth's equator. 
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2. The systematic error of gravitational origin in the LAGEOS-LAGEOS II 

test 



The reahstic eva^ 



Ciu: 



2006 



'olini and Pav 



lorij 



2005a 



IS 



uation of the total error bud get of the LAGEOS- 



2006a 



2004) raised a lively debate 



2007c 



Lucchesi 



(ICiufolini and Pavlis 



AGEOS II node test 



2005 



Ciufolini et al 



20051 ). mainly focussed on the impact of the 



static and time-varying parts of the Newtonian component of the Earth's gravitational 
potential through the aliasing secular precessions induced on a satellite's node. 

In a realistic scenario the path of a probe is not only affected by the gravitomagentic 
field but also by a huge number of other competing orbital perturbations of gra vitational and 



non-g ravitational origin. The most important non-conservative accelerations (IMilani et al 



19871 ) are the direct solar radiation pressure, the Earth's albedo and various subtle thermal 



effects depending on the the physical 



state ( 



Inversi and Vespe 



Lucchesi et al 



2004 



1994 



Vespe 



Ries et al 



properties of the sate l 



1999 



Lucchesi 



2001 



ite's surface and its rotational 



Lucchesi 



2002 



2003 



2004 



2003al ): however, the nodes of LAGEOS-like satellites are 
sensitive to them at a ^ 1% level only. Much more important is the impact that the 
oblateness of the Earth, due to its diurnal rotation, has on the satellite's dynamics. Indeed, 
the most insidious perturbations are those induced by the static part of the Newtonian 



component of the multipolar expansion in spherical harmonics 



Jf 



2,4,6,... of the 



19661 ): indeed, they affect the 



gravitational potential of the central rotating mass (iKaulal 
node with effects having the same signature of the relativistic signal of interest, i.e. linear 
trends which are orders of magnitude larger and cannot be removed from the time series of 
data without affecting the Lense-Thirring pattern itself as well. The only thing that can 
be done is to model such a corrupting effect as most accurately as possible and assessing. 



^The relation among the even zonals Jg and the normalized gravity coefficients Ceo which 
are customarily determined in the Earth's gravity models, is Jg = —y/2i+ 1 Ceo. 
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reliably and realistically, the impact of the residual mismodelling on the measurement 
of the frame-dragging effect. The secular precessions induced by the even zonals of the 
geopotential can be written as 

Qseopot ^J2Q,eJ,, (5) 

where the coefficients = 2,4,6, ... depend on the parameters of the Earth {GM and 
the equatorial radius R) and on the semimajor axis a, the inclination i and the eccentricity 
e of the satellite. For example, for £ = 2 we have 



3 fR 
'-n — 
2 V a 



cosz 



:i-e 



2^2 ' 



(6) 



n 



^GUja? is the Keple rian mean motion. They have been analytically computed up to 



20 in, e.g., (jlorio 



20031 ). Their mismodelling can be written as 



^j^geopot <YJ^.i bJl, (7) 
1=2 

where 8Ji represents our uncertainty in the knowledge of the even zonals Jn 

A common feature of all the competing evaluations so far published is that the 
systematic bias due to the static component of the geopotential was calculated always by 
using the released (more or less accurately calibrated) sigmas oj^ of one Earth gravity 
model solution at a time for the uncertainties bJn in the even zonal harmonics, so to say 
that the model X yields a x% error, the model Y yields a y% error, and so on. 

Since a trustable calibration of the formal, statistical uncertainties in the estimated 
zonals of the covariance matrix of a global solution is always a difficult task to be 
implemented in a reliable way, a much more realistic and conservative approach consists, 
instead, of taking the difference! AJ^ of the estimated even zonals for different pairs of 
Earth gravity field solutions as representative of the real uncertainty 63^ in the zonals 



^See Fig. 5 of (ILucchesi 



20071 ) for a comparison of the estimated C40 in different models. 



- 13 - 



( iLerch et al I Il994l ). In Table [Tl-Table [TO] we present our results for the most recent 
GRACE-based models released so far by different institutions and retrievable on the 
Internet at'^ 



http: / /ic gem.gfz-potsdam.de /ICGEM /ICGEM.htm: 



EIGEN-GRACE02S (IReigber et al 



GFZ (Potsdam, Germany), GGM02S flTaplev et al 



2005h and EIGEN-C G03C flForste et al 



The mode 



2nn5h and GGM 03S dTaplev et al 



s used are 



2005) from 



from CSR (Austin , Texas), ITG-Grace02s (jMayer-Giirr et al I l2006l ) and ITG-Grace03s 



20071 ) 



( iMayer-Giirr 



2007h from IGG fBo nn. Germany), JEM01-RL03B from JPL (NASA, USA) 



and EGM2008 ( 
followed also by 



'avlis et a. 



20081) from NGA (USA). Note that this approach was explicitly 



Ghifohnil (jl996h with the JGM3 and GEMT-2 models. 



The systematic bias evaluated with a more realistic approach is about 3 to 4 times 



larger than one can obtain by only using this or that particular mode 



The scatter is still 



20041 ). In particular. 



quite large and far from the 5 — 10% claimed in (jCiufolini and PavlisI 
it appears that J4, Jg, and to a lesser extent Jg, which are just the most relevant zonals 
for us because of their impact on the combination of eq. ([3]), are the most uncertain 
ones, with discrepancies AJ^ between different models, in general, larger than the sum 
of their sigmas cij^, calibrated or not. This is an important feature because the other 



alternative combinations p roposed involving more satellites (jlorio and Doornbos 



2005 



Vespe and Rutigliano 



2OO5I ) should be less affected since they cancel out the impact of J4 



and Je as well. 

Another approach that could be followed to take into account the scatter among the 
various solutions consists in computing mean and standard deviation of the entire set of 
values of the even zonals for the models considered so far, degree by degree, and taking 
the standard deviations as representative of the uncertainties 5Ji,^ = 4,6,8, .... It yields 



I thank J Rie s, CSR, and M Wa tkins (JPL) for having provided me with the even zonals 



of the GGM03S (iTaplev et al 



20071) and JEM01-RL03B models. 
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5^l = 15%. 



It must be recalled that also the further bias due to the cr oss-coupling; between J2 and 



the orbit inclination, evaluated to be about 9% in (Iloridl2007d ). must be added. 



3. A conservative evaluation of the impact of the geopotential on the LARES 

mission 

The combination which will be us ed for meas uring the Lense-Thirring effect with 



LAGEOS, LAGEOS II and LARES is fllorio 



2005ah 



^LAGEOS ^ ^^^LAGEOS II ^ ^^J^LARES^ 



The coefficients ki and k2 entering eq. (IHl) are defined as 

oLARESoLAGEOS oLAGEOSoLARES 

U ".2 ".4 l\2 1^4 n Qcqy 

"'I qLAGEOS IlQLARES„f2LARESQLAGEOS II U.OU J ( , 



oLAGEOSoLAGEOS II oLAGEOS IIoLAGEOS 

u ii2 Ha ii2 lii fi fifilQ 

"^2 qLAGEOS IIfjLARES_f2LARESfjLAGEOS II U.UUIJ. 



(9) 



The combination eq. ([8]) cancels out, by construction, the impact of the first two even 
zonals. The total Lense-Thirring effect, according to eq. ([8]) and eq. iQ, amounts to 
47.8 mas yr~^. The systematic error due to the uncancelled even zonals Jq, Jg, ... can be 
conservatively evaluated as 

5/i < 5^ |f2L/GEOS ^ ^^^LAGEOS II ^ ^^l^L/RES | ^j^ (^g) 
£=6 

Of crucial importance is how to assess 5Je. By proceeding as in Section [2] and by using 
the same models up to degree i = 60 because of the lower altitude of LARES with respect 
to LAGEOS and LAGEOS II which brings into play more even zonals, we have the results 
presented in Table [TTl 
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Table 1: Impact of the mismodelling in the even zonal harmonics on = 
j^LAGEOs ^ ^^j^LAGEOs II ^j^^ £ = 4, . . . , 20, in mas yr-i. Recall that Ji = -VWTT Cm] 



for the uncertainty in the even zonals we h ave taken here the d ifference AC, 



between the model X = E 



GRACE02S (iReigber et al 



GEN- CG03C f lForste et al 



C 



(x) 



C 



(Y) 



20051) and the model Y = EIGEN- 



20051). EIGEN-CG03C combines data from CHAMP (860 days 
out of October 2000 to June 2003), GRACE (376 days out of February to May 2003, 
July to December 2003 and February to July 2004) and terrestrial measurements; EIGEN- 
GRACE02S is based on 110 days (out of August and November 2002 and April, May and 
August 2003) of GRACE-only CPS-GRACE high-low satellite-to-satellite data, on-board 
measurements of non-gravitational accelerations, and especially GRACE intersatellite track- 
ing data. (Tx/Y are the covariance calibrated errors for both models. Values of fi smaller 
than 0.1 mas yr~^ have not been quoted. The Lense-Thirring precession of the combination 
of eq. ([3]) amounts to 47.8 mas yr~^. The percent bias 5/i has been computed by normal- 
izing the linear sum of Z^, £ = 4, . . . , 20 to the Lense-Thirring precession. The discrepancies 
between the models are significant since AC^o are larger than the linearly added sigmas for 





4,. ..16. 








ACm (EIGEN-CG03C-EIGEN-GRACE02S) 


CTX + cry 


fi (mas yr-i) 


4 


1.96 X 10^11 


1.01 X 10-11 


7.3 


6 


2.50 X 10-11 


4.8 X 10-12 


5.4 


8 


4.9 X 10-12 


3.3 X 10-12 


0.2 


10 


3.7 X 10-12 


3.4 X 10-12 




12 


2.5 X 10-12 


2.3 X 10-12 




14 


6.1 X 10-12 


2.1 X 10-12 




16 


2.1 X 10-12 


1.7 X 10-12 




18 


6 X 10-13 


1.7 X 10-12 




20 


1.7 X 10-12 


1.7 X 10-12 




total bias = 27% 
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able 2: Impact of t he mismodelling in the e ven zonal harmonics as so lved for in X=GGM02S 



( jTapley et al 



20051) and Y=ITG-Grace02s f Maver-Giirr et al 



2006|). GGM02S is based on 



363 days of GRACE-only data (GPS and intersatellite tracking, neither constraints nor 
regularization applied) spread between April 4, 2002 and Dec 31, 2003. The a are formal 
for both models. AC^o are always larger than the linearly added sigmas, apart from £ = 12 



and 


£ = 18 










£ 




(GGM02S-ITG-Grace02s) 


CTX + CTy 




fn (mas yr-i) 


4 


1.9 X 


10-" 


8.7 X 10- 


-12 


7.2 


6 


2.1 X 


10-" 


4.6 X 10- 


-12 


4.6 


8 


5.7 X 


10-12 


2.8 X 10- 


-12 


0.2 


10 


4.5 X 


10-12 


2.0 X 10" 


-12 




12 


1.5 X 


10-12 


1.8 X 10- 


-12 




14 


6.6 X 


10-12 


1.6 X 10- 


-12 




16 


2.9 X 


10-12 


1.6 X 10" 


-12 




18 


1.4 X 


10-12 


1.6 X 10- 


-12 




20 


2.0 X 


10-12 


1.6 X 10- 


-12 





total bias 5/i = 25% 
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able 3: Impact of t 



( jTapley et al 



le mismodelling in the even zonal harmonics as solved for in X=GGM02S 



20051 ) and Y=EIGEN-CG03C flF5rste et al 1120051 ). The a are formal for 
GGM02S, calibrated for EIGEN-CG03C. AC^o are always larger than the linearly added 
sigmas. 





AC^o (GGM02S-EIGEN-CG03C) 


(Tx + CXy 


(mas yr-i) 


4 


1.81 X 10^11 


3.7 X 10-12 


6.7 


6 


1.53 X 10"^^ 


1.8 X 10-12 


3.3 


8 


1.5 X 10^12 


1.1 X 10-12 




10 


4.9 X 10-^2 


8 X 10-13 




12 


8 X 10-13 


7 X 10-13 




14 


7.7 X 10-12 


6 X 10-13 




16 


3.8 X 10-12 


5 X 10-13 




18 


2.1 X 10-12 


5 X 10-13 




20 


2.3 X 10-12 


4 X 10-13 





total bias = 22% 
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able 4: Bias due to the mismodelling in the even zonals of the models X=ITG-Grace03s 



( iMayer-Giirr 



20071 ). based on GRACE-only accumulated normal equations from data out 



of September 



2002-April 2007 (n either apriori information nor regularization used), and 



Y=GGM02S (iTaplev et al 



20051 ). The a for both models are formal. AC^o are always 



larg 


er than the linearly added sigmas, apart from (. 


= 12 and i = 18 


£ 


ACeo (ITG-Grace03s-GGM02S) 




fi (mas yr-i) 


4 


2.58 X 10-11 


8.6 X 10-12 


9.6 


6 


1.39 X 10-11 


4.7 X 10-12 


3.1 


8 


5.6 X 10-12 


2.9 X 10-12 


0.2 


10 


1.03 X 10-11 


2 X 10-12 




12 


7 X 10-13 


1.8 X 10-12 




14 


7.3 X 10-12 


1.6 X 10-12 




16 


2.6 X 10-12 


1.6 X 10-12 




18 


8 X 10-13 


1.6 X 10-12 




20 


2.4 X 10-12 


1.6 X 10-12 




total bias 5/i = 27% 
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able 5: Bias due to the mismodelling in the even zonals of the models X = GGM02S 



( jTapley et al 



20051 ) and Y = GGM03S (jTapley et al 1 120071 ) retrieved from data spanning 



January 2003 to December 2006. The cr for GGM03S are calibrated. AC^o are larger than 
the linearly added sigmas for £ = 4,6. (The other zonals are of no concern) 





ACa (GGM02S-GGM03S) 


cTx + o-y 


fi (mas yr-i) 


4 


1.87 X IQ-^^ 


1.25 X 10-11 


6.9 


6 


1.96 X 10-11 


6.7 X 10-12 


4.2 


8 


3.8 X 10^12 


4.3 X 10-12 


0.1 


10 


8.9 X 10-12 


2.8 X 10-12 


0.1 


12 


6 X 10-13 


2.4 X 10-12 




14 


6.6 X 10-12 


2.1 X 10-12 




16 


2.1 X 10-12 


2.0 X 10-12 




18 


1.8 X 10-12 


2.0 X 10-12 




20 


2.2 X 10-12 


1.9 X 10-12 




total bias 5ix = 24% 
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Table 6: Bia s due to the mismo delling; in the even zonals of the mod els X = EIGEN 



GRACE02S flReigber et al 



2005|) and Y = GGM03S (ITaplev et al 1120071 ). The a for both 



models are calibrated. AC^o are always larger than the linearly added sigmas apart from 






14,18. 








AC^o (EIGEN-GRACE02S-GGM03S) 


(Tx + o-y 


(mas yr-i) 


4 


2.00 X 10^11 


8.1 X 10-12 


7.4 


6 


2.92 X 10^1^ 


4.3 X 10-12 


6.3 


8 


1.05 X 10^11 


3.0 X 10-12 


0.4 


10 


7.8 X 10-12 


2.9 X 10-12 


0.1 


12 


3.9 X 10-12 


1.8 X 10-12 




14 


5 X 10-13 


1.7 X 10-12 




16 


1.7 X 10-12 


1.4 X 10-12 




18 


2 X 10-13 


1.4 X 10-12 




20 


2.5 X 10-12 


1.4 X 10-12 





total bias = 30% 



- 21 - 



Table 7: Bias due to the mismodelling in the even zonals of the models X = J 



based on 49 months of GRACE-only data, and Y = GGM03S flTaplev et al 



;M01-RL03B, 



2007h . The a 



for GGM03S are calibrated. AC^o are always larger than the linearly added sigmas apart 



from 


£=16. 






^ 


AC^o (JEM01-RL03B-GGM03S) 


CTx + cry 


fi (mas yr-i) 


4 


1.97 X 10-11 


4.3 X 10-12 


7.3 


6 


2.7 X 10^12 


2.3 X 10-12 


0.6 


8 


1.7 X 10^12 


1.6 X 10-12 




10 


2.3 X 10-12 


8 X 10-13 




12 


7 X 10-13 


7 X 10-13 




14 


1.0 X 10-12 


6 X 10-13 




16 


2 X 10-13 


5 X 10-13 




18 


7 X 10-13 


5 X 10-13 




20 


5 X 10-13 


4 X 10-13 





total bias 5/1 = 17% 
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Table 8: Bias due to the mismodelline; in t 



le even zonals of the models X = JEM01-RL03B 



and Y = ITG-Grace03s ( Mayer-Gun 



20071 ). The a for ITG-Grace03s are formal. AC^o are 



always larger than the linearly added sigmas. 





AC^o (JEM01-RL03B-ITG-Grace03s) 


c^x + cfy 


ft. (mas yr-i) 


4 


2.68 X 10-11 


4 X 


10-13 


9.9 


6 


3.0 X 10-12 


2 X 


10-13 


0.6 


8 


3.4 X 10-12 


1 X 


10-13 


0.1 


10 


3.6 X 10-12 


1 X 


10-13 




12 


6 X 10-13 


9 X 


10-14 




14 


1.7 X 10-12 


9 X 


10-14 




16 


4 X 10-13 


8 X 


10-14 




18 


4 X 10-13 


8 X 


10-14 




20 


7 X 10-13 


8 X 


10-14 





total bias 5/x = 22% 
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Table 9: Aliasing effect of the mism odelling; in the even zonal harmo nics estimated in the 



X=ITG-Grace03s flMaver-Giirr 



20071 ) and the Y=EIGEN-GRACE02S flReigber et al 1120051 ) 
models. The covariance matrix a for ITG-Grace03s are formal, while the ones of EIGEN- 
GRACE02S are calibrated. AC^o are larger than the linearly added sigmas for 
apart from £ = 18. 



4,. ..,20, 





KCm (ITG-Grace03s-EIGEN-GRACE02S) 


o'x + cry 


fi (mas yr-i) 


4 


2.72 X 10^11 


3.9 X 10-12 


10.1 


6 


2.35 X 10^11 


2.0 X 10-12 


5.1 


8 


1.23 X 10^11 


1.5 X 10-12 


0.4 


10 


9.2 X 10-12 


2.1 X 10-12 


0.1 


12 


4.1 X 10-12 


1.2 X 10-12 




14 


5.8 X 10~i2 


1.2 X 10-12 




16 


3.4 X 10-12 


9 X 10-13 




18 


5 X 10-13 


1.0 X 10-12 




20 


1.8 X 10-12 


1.1 X 10-12 





total bias 5/i = 37% 
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Table 10: Im pact of the mism odelling: in the even zonal har monics estimated in the 



X=EGM2008 f Pavhs et al. 



120081) and the Y=EIGEN-GRACE02S flReigber et al 1120051 ) mod- 
els. The covariance matrix a are calibrated for both EGM2008 and EIGEN-GRACE02S. 
ACfo are larger than the linearly added sigmas for £ = 4, 20, apart from I = 18. 



i 


ACio (EGM2008-EIGEN-GRACE02S) 


(Tx + (Ty 




fi (mas yr-i) 


4 


2.71 X 10^11 


8.3 X 10- 


-12 


10.0 


6 


2.35 X 10"^^ 


4.1 X 10- 


'12 


5.0 


8 


1.23 X 10-11 


2.7 X 10- 


-12 


0.4 


10 


9.2 X 10-12 


2.9 X 10- 


-12 


0.1 


12 


4.1 X 10-12 


1.9 X 10- 


-12 




14 


5.8 X 10-12 


1.8 X 10- 


-12 




16 


3.4 X 10-12 


1.5 X 10- 


-12 




18 


5 X 10-13 


1.5 X 10- 


-12 




20 


1.8 X 10-12 


1.5 X 10- 


-12 





total bias Sfj, = 33% 
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Table 11: Systematic percent error 6fi in the measurement of the Lense-Thirring effect with 
LAGEOS, LAGEOS II and LARES according to eq. and 6Je = AJe up to degree i = QO 
for the global Earth's gravity solutions considered here. 



Models compared 

EGM2008-JEM01-RL03B 8% 

EGM2008-GGM02S 27% 

EGM2008-GGM03S 5% 

EGM2008-ITG-Grace02 4% 

EGM2008-ITG-Grace03 0.1% 

EGM2008-EIGEN-CG03C 38% 

EGM2008-EIGEN-GRACE02S 53% 

JEM01-RL03B-GGM02S 28% 

JEM01-RL03B-GGM03S 10% 

JEM01-RL03B-ITG-Grace02 8% 

JEM01-RL03B-ITG-Grace03s 8% 

JEM01-RL03B-EIGEN-CG03C 44% 
JEM01-RL03B-EIGEN-GRACE02S 57% 

GGM02S-GGM03S 24% 

GGM02S-ITG-Gracc02 28% 

GGM02S-ITG-Gracc03s 26% 

GGM02S-EIGEN-CG03C 27% 

GGM02S-EIGEN-GRACE02S 36% 

GGM03S-ITG-Grace02 5% 

GGM03S-ITG-Grace03s 5% 

GGM03S-EIGEN-CG03C 36% 

GGM03S-EIGEN-GRACE02S 52% 

ITG-Grace02-ITG-Graee03s 4% 

ITG-Graee02-EIGEN-CG03C 39% 

ITG-Graee02-EIGEN-GRACE02S 54% 

ITG-Grace03s-EIGEN-CG03C 38% 
ITG-Grace03s-EIGEN-GRACE02S 53% 
EIGEN-CG03C-EIGEN-GRACE02S 27% 
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It must be stressed that they may be still optimistic: indeed, computations for i > QQ 
become unreliable because of numerical instability of the results (obtained with two different 
softwares) . 

If, instead, one assumes 6Je = S£, i = 2,4,6, ... i.e., the standard deviations of the sets 
of all the best estimates of Je for the models considered here the systematic bias, up to 
i = 60, amounts to 26%. Again, also this result may turn out to be optimistic for the same 
reasons as before. 

It must be pointed out that the evaluations presented here rely upon cal culations of th e 



Kaulal fll966h 



coefficients performed with the well known standard approach by Kaula 
it would be important to try to follow also different computational strategies in order to 
test them. 

It is worthwhile noting that also the impact of the subtle non-gravitational perturbations 
will be different with respect to the original proposal because LARES will fly in a different 
and lower orbit and its thermal behavior will probably be different with respect to LAGEOS 
and LAGEOS II. The reduction of the impact of the thermal accelerations, like the 
Yarkowsky-Sch ach effec t s, sho uld have been reached with two concentric spheres. However, 



as explained in lAndresI (j2007l ). this solution will increase the floating potential of LARES 
because of the much higher electrical resistivity and, thus, the perturbative effects produced 
by the charged particle drag. Moreover, drag will increase also because of the lower orbit of 
the satellite, both in its neutral and charged components. Also the Earth's albedo, with its 
anisotropic components, should have a major effect. 

Another point which must be considered is the realistic orbit accuracy obtainable for 
LARES. Indeed, at a lower orbit the normal points RMS will be probably higher with 
respect to the present RMS obtained for the two LAGEOS satellites (a few mm), as we 
presently know for the Stella and Starlette normal points. Of course, such an accuracy is a 
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function of several aspects. 

4. Conclusions 

In this Chapter we have shown how the so far published evaluations of the total 
systematic error in the Lense-Thirring measurement with the combined nodes of the SLR 
LAGEOS and LAGEOS II satellites due to the classical node precessions induced by 
the even zonal harmonics of the geopotential are optimistic. Indeed, they are all based 
on the use of the covariance matrix's sigmas, more or less reliably calibrated, of various 
Earth gravity model solutions used one at a time separately in such a way that the model 
X yields an error of x%, the model Y yields an error y%, etc. Instead, comparing the 
estimated values of the even zonals for different pairs of models allows for a much more 
realistic evaluation of the real uncertainties in our knowledge of the static part of the 
geopotential. As a consequence, the bias in the Lense-Thirring effect measurement is about 
3 — 4 times larger than that so far claimed, amounting to various tens percent (37% for the 
pair EIGEN-GRACE02S and ITG-GRACE03s, about 25 - 30% for the other most recent 
GRACE-based solutions). 

Applying the same strategy to the ongoing LARES mission shows that the goal of 
reaching a 1% measurement of the Lense-Thirring effect with LAGEOS, LAGEOS II and 
LARES is optimistic. Indeed, since LARES will orbit at a lower altitude with respect to the 
LAGEOS satellites, a larger number of even zonal harmonics are to be taken into account. 
Assessing realistically their impact is not easy. Straightforward calculations up to degree 
i = 60 with the standard Kaula's approach yield errors as large as some tens percent. Such 
an important point certainly deserves great attention. Another issue which may potentially 
degrade the expected accuracy is the impact of some non-gravitational perturbations which 
would have a larger effect on LARES than expected because of its lower orbit. 
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